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Introduction 


In Unit 14 on vector algebra we showed how a physical quantity could be 
categorized as either a scalar or a vector, and we spent some time on the algebra 
of vectors, i.e. adding, scaling and multiplying vectors. A scalar quantity has 
magnitude only and can be represented by a real number; for instance, mass and 
temperature are scalars. A vector quantity has two properties associated with it, its 
magnitude and a definite direction in space. We represent a vector on a diagram 
by an arrow or algebraically by components. Velocity and force are examples of 
physical quantities which can be represented by vectors. 


In Unit 14 the values of the scalars and the components of the vectors did not 
depend on where they were measured. 


In this unit we shall consider scalars and vectors which have different values at 
different points in a region of space. For instance, you know from practical 
experience that the temperature in a bowl of hot soup varies quite significantly 
from the middle to the edge of the bowl. That is why we try to drink the soup 
from the edge rather than from the middle where we are likely to get a burnt 
mouth! Temperature is a scalar quantity and its value depends on where in the 
soup it is measured. The temperature distribution is an example of a scalar field. 


Similarly, the force on a nail in the region surrounding a magnet will depend on 
how near to the magnet we bring the nail. The magnetic force field of the magnet 
has a definite magnitude and direction which varies with the position of the point 
at which it is measured. This is an example of a vector field. 


This unit is about scalar and vector fields. 


Study guide 


This unit is divided into four sections and they should be studied in the order in 
which they appear. The audio-tape is part of Subsection 3.3 and the television 
programme should be viewed at the beginning of Section 4. If you have not had 
time to study the first three sections before watching the television programme 
then a quick skim through Sections | and 3 should be sufficient for you to follow 
and enjoy the programme; then return and study the unit in more detail. 


In Section 1 we define scalar and vector fields and show how we can describe 
them geometrically. In Section 2 we define the spatial rate of change of a scalar 
field called the gradient of that scalar field. Section 3 introduces a new kind of 
integral of a vector field called the line integral. You should listen to the tape for 
this unit while studying this section; it takes you through several line integrals to 
help you gain confidence at evaluating them. We then, in Section 4, define a new 
vector field related to the spatial rate of change of the old vector field and called 
its curl. This definition uses the line integral introduced in Section 3. The television 
programme, which is associated with this section, shows you the physical 
significance of curl by considering the velocity vector fields of three fluid flows. 


Finally, a note to those students who have previously studied MS283. You should 
be familiar with the material contained in this unit. You are recommended to 
work through the exercises and to revise those sections whose exercises you find 
difficult. 
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1 Scalar and vector fields 


1.1 A scalar field defined 


Figure | shows the temperature at different points of a flat plate when unevenly 
heated and cooled. The value of the temperature depends on where on the plate it 
is measured. If we define Cartesian coordinate axes Ox and Oy as shown, then the 
variable T representing the temperature will be a function of x and y; ie, 

T= T(x, y). This function of two variables is just an example of the concept of a 
function of more than one variable introduced in Unit 25. Figure 1 


On 
heat source 


As another example, Figure 2 shows the (atmospheric) pressure distribution over 
the U.K. at midday on 30 September 1980. 


Figure 2 


The solid lines are called isobars and join up places where the atmospheric 

pressure is the same. Atmospheric pressure at a given time is a scalar quantity and 

its value depends on where in the U.K. it is measured. For instance, at midday on 

30 September 1980, in Manchester it was 1020 millibars and in Penzance it was 

1028 millibars. 1 millibar = 100 Nm~? 


If we choose London as an origin and set up a Cartesian coordinate system with 
Ox pointing East and Oy pointing North, then each place in the U.K. can be 
represented by its coordinates (x,y). Then the variable P representing the 
atmospheric pressure is a function of two variables x and y; ie, P = P(x,y). 


These two examples are pictorial representations of scalar fields. In each of the 
examples above we have: (i) a scalar quantity (e.g. temperature or atmospheric 
Pressure), (ii) a region of space over which each scalar is defined (e.g. the flat plate 
or the U.K.), and (iii) the values of the scalar quantity at different points of the 
region. 


In each of the above cases, the variable representing the scalar quantity is a 
function of the two variables x and y; these are two-dimensional scalar fields. The 
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temperature distribution in the room where you are reading this unit would be an 
example of a three-dimensional scalar field, which we now define: 


If @ is a function which associates a unique scalar with each point Waican sefinon two. 3 


in a given region of space, then @ is called a scalar field function similarly, but then the domain 
or simply a scalar field, is a region of the plane and 
not a region of space. 


The domain of the function ¢ is the given region of space and the codomain is the 
set of all scalars. 


Example 1 


The light intensity from a point source of light, such as a small light bulb, is 
inversely proportional to the square of its distance from the bulb. The intensity at 
distance 1 metre is Jy. Specify a scalar field function which models the light 
intensity and state the region of space over which it is defined. 


Solution 
Suppose that we represent the light source by a point at the origin, and the light 
intensity at a point distance r metres from the light source by J; then I a = and 


Timp when = ip Thus =e? 


iets 


The scalar field function is J = 3 ( in Cartesian coordinates). The 


To 
ae +y? +27) 
function is not defined when r = 0 and so the region on which the scalar field is 
defined is the whole of space with the origin (0,0,0) removed. Note that this is a y 
three-dimensional scalar field. 


Exercise 1 
A thin circular sheet of metal of radius 4 metres has a concentric hole of radius 1 metre cut (A) 
out of it (see Figure 3). The inner perimeter is maintained at a temperature of 0°C and the x 


outer perimeter is maintained at a temperature of 100°C. At a point P on the sheet, given 
by position vector r, the temperature is proportional to log, r, where r = |r|. Specify the 
scalar field function which models the temperature and state the region of the plane over 
which it is defined, 


[Solution on p. 44] 


Figure 3 


1.2 Contour curves and contour surfaces 


The scalar fields in the last subsection were presented in three different ways. For 
the first scalar field of temperature over a plate, as given in Figure 1, we 
represented the data as values of the temperature at a few points on the plate. For 
the second scalar field of atmospheric pressure in the U.K., as given in Figure 2, 
we represented the data by isobars (or lines joining points at which the 
atmospheric pressure is the same) on a map of the U.K. (In both cases this is not 
a complete representation of the scalar field because for instance, there is a 
temperature associated with every point of the flat plate, and not just at the points 
shown in Figure 1.) Finally, in Example 1 and Exercise 1 the scalar fields were 
given by a formula of some type. 


Of these modes of representation, the isobar example gives a very good way of 
visualizing pictorially the physical situation. Figure 4 shows the temperature 
distribution over an unevenly heated flat plate similar to that shown in Figure 1. 
However instead of giving the value of the temperature at different points as we 
did before, we have drawn curves through those points on the plate where the 
temperature is the same. Of course we cannot show all possible curves because 
that would require a curve through every point on the plate, so we choose certain 
values of the temperature and the curves give us a pictorial idea of the scalar field. _heat source 
So, for example, for each point on the curve AB the temperature equals 300°C. Figure 4 
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For a temperature field the curves are called isotherms, and in this example we 
have drawn the isotherms for multiples of 50°C. Notice that in one part of the 
plate the isotherms are much closer together than for the rest of the plate; the 
temperature is changing more rapidly in the region where the isotherms are close 
together. 


We now extend this pictorial representation for pressure and temperature to more 
general two-dimensional scalar fields. Suppose we have a two-dimensional scalar 
field (x, y) defined over a region R. We can represent this scalar field pictorially 
by drawing curves through those points of R which have equal values of . Figure 
5, for example, shows four such curves where C,, C,, C3 and C, are constants. 
The curves that we have drawn are called contour curves. 


For example, consider the scalar field 

g=x+y? 
defined over the rectangle —1 < x < 1, —2 < y <2. The contour curves are 
defined by 

@=constant, or x? + y*= constant. 


These curves are circles or parts of circles drawn within the rectangular region for 
which @ is defined, Figure 6 shows the contour curves 
(x,y) = 0.25; (x,y) = 1; (x,y) = 2.25; (x,y) = 4. 


Again we cannot draw all possible contour curves, only some of them, 


Exercise 2 
Draw some of the contour curves to represent pictorially the scalar field ¢ = xy defined 
over the disc x? + y? <4. 


[Solution on p. 44) 


You are probably familiar with the word ‘contour’ from Ordnance Survey maps. 
These are lines on a map joining up places which are at the same height above sea 
level. The contours can give us a lot of information about the lie of the land, Here 
is an example; see if you can answer Exercise 3, just from the contour curves. 


Figure 7 


Exercise 3 
(i) Where would you expect to find (a) a hill, and (b) a lake? 


(ii) Where do you think that the land is fairly level? 


(iii) Suppose you set off from point A and walked in the direction of the arrow, would you 
expect an easy walk or a hard climb? 


[Solution on p. 44] 


This exercise illustrates one important feature of a scalar field which can be 
obtained from its contour curves. The rate of increase (or decrease) of a scalar 
field can be deduced from the closeness of the contour curves. If the contour 
curves are close together then the scalar field is changing quite rapidly in that part 
of the region. 


Figure 6 
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So far in this subsection we have considered two-dimensional fields. Now we turn 
our attention to three-dimensional scalar fields. A pictorial representation of a 
two-dimensional scalar field #(x, y) is given by the contour curves defined by z 
(x, y) = constant. We can represent a three-dimensional scalar field (x, y,z) in an | 
analogous way. If we join up points for which the three-dimensional scalar field 

(x,y,z) has the same value, we obtain the surface with equation $(x, y,z) = C 
where C is constant. The family of surfaces given by (x,y,z) = C, for different 
values of the constant C, are called the contour surfaces of the scalar field and they 
give a pictorial representation of the scalar field. 


For example, the contour surfaces of the scalar field J = Io/(x? + y? + 2?) given in 
Example 1 are concentric spheres with centres at the origin. Two of the contour 
surfaces, for J = 4 and J = 9, are shown in Figure 8. 


1.3 A vector field defined 


Figure 9 shows what happens to iron filings when placed near a bar magnet. The 
iron filings have aligned themselves in a symmetrical pattern and they are 
experiencing a force due to magnetism. If we place small compasses at various 
points near the magnet we find the arrow on a compass orientates itself along the 
iron filings. The compasses show us the direction of the magnetic field at each 
point. The magnitude of the field decreases as we go further away from the 


magnet. The arrows in Figure 10 represent the magnetic field due to the bar \ 
magnet at various points near to the magnet. At each point in a region close to Pee 
the magnet there is a unique vector called the magnetic field at that point. This is el: 
an example of a vector field. Ses © eA: 

. 
Here is another example. The arrows in Figure 11 represent the surface velocity at i v 
various points in a river. eaad 

- 


Figure 10 


Figure 11 


Near the river bank the water is almost at rest, so the arrows are somewhat 
shorter in length (indicating a smaller velocity magnitude near the edge of the 
river). The velocity at a particular point is the velocity that a small floating object, 
such as a leaf, would have at that point. At each point on the river’s surface there 
is a unique velocity; this is an example of a velocity vector field. 


You have seen something like this before in Unit 2 where you met direction fields. 
These were a distribution of orientations in space, i.e. directions without a sense 
being attached, but no magnitudes were associated with the elements of the 
direction field. We now have a distribution of vectors in a region of space, i.e. 
quantities with magnitude and direction, which may take on different values at each 
point in the region. This leads to the following definition: 


If F is a function which associates a unique vector with every point 
in a given region of space, then F is called a vector field function 
or simply a vector field, 
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The domain of the function is the given region of space and the codomain is the 
set of all vectors. To describe a vector field completely we require the following: 

(i) the region of space on which the vector field is defined, and (ii) the magnitude 
and direction of the vector at each point of the region. 


Example 2 


The Earth’s gravitational field can be considered as an example of a vector field: 
as you will see in Unit 30, the force on a body which is at a distance r from the 
centre of the Earth has magnitude inversely proportional to the square of the 
distance r and is directed towards the Earth’s centre. Specify the vector field and 
the region of space over which it is defined. 


Solution 
The vector quantity is the gravitational force acting on the body due to the Earth. 
This force acts at points (x, y,z) in space such that 


Je+y+272R 
where R is the radius of the Earth, and x, y, z are Cartesian coordinates with 
origin at the centre of the Earth. The vector field function is — Fee, where e, is a 
unit vector whose direction is from O to P (see Figure 12), The magnitude F, is 
proportional to 1/(OP)?, so 
pee c 
* OP? (x? +? + 27) 
for some constant C, and the direction of F is along the vector 
PO = — (xi + yj + zk) 
so that 
mae 
oP =€ (xi + yj + 2k) 
=k “@+ypte) joypad 
[oP GP +P +2) (ey yee 
_ —C(xi + yj + 2k) 
= @+yPtepe: 


The region over which the vector field F is defined is the set of points (x, y, z 
outside the sphere x? + y? + 2? = R?, 

This is an example of a three-dimensional vector field because the function F is 
defined on a three-dimensional region of space. 


The components of the vector field function relative to a given Cartesian 
coordinate system are themselves scalar field functions defined over the same 
region of space as F, Thus if F is a three-dimensional vector field function defined 
on a region D and 

F(x, y,2) = F(x, y,2)i + Fa(x,y,2)) + Fy(x. 9. 2)k 
then F (x,y,z), Fa(x,y,z) and F(x, y,2) are scalar field functions also defined on D. 


The following exercise gives an example of a two-dimensional vector field. 


Exercise 4 


A flat horizontal metal circular plate of radius 10cm has a small hole in the middle through 
which is injected water at a constant rate. It is found that the velocity of the water at any 
point P on the plate is directed outwards along the line joining P to the centre of the plate 
O and has magnitude inversely proportional to the distance from 0. Specify the vector field 
and the region over which it is defined, in terms of Cartesian coordinates with origin at 0. 


{Solution on p. 44] 


1.4 Vector field lines 


The arrows on the diagrams in the last subsection are one way of describing a 
vector field pictorially. The pictures are however incomplete because there is a 


Figure 12 
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vector quantity associated with every point in the region. A more complete picture 
can be obtained by drawing continuous curves in the region such that at any point 
the tangent to the curve is parallel to the direction of the vector field function at 
that point. These curves are called the (vector) field lines of a vector field. For 
instance, the orientation of iron filings near a bar magnet suggests a family of 
continuous curves. Two such curves are shown in Figure 13. At any point on one 
of the curves the orientation of a compass needle would be along the tangent to 
the curve; one is shown in the figure. 


Exercise 5 
Sketch the field lines for the Earth’s gravitational force field. 
[Solution on p. 44) 


When the vector field represents a force the field lines are often called lines of 
force. 


Here is another example of field lines. The arrows in Figure 14 represent the 
velocity of the surface ofa river at various points. If we drop some leaves onto the 
river they will float. If we draw a line showing the path of one of these leaves, then 
at each point on the path, the tangent to the path is parallel to the direction of the 
velocity vector. This line is a field line for the velocity vector field, For the motion 
of a liquid the field lines are often called flow lines. 


You have met field lines before in Unit 2 of this course. The trajectories of a 
direction field are analogous to the field lines of a vector field. 

The field lines give a convenient way of describing a vector field pictorially, and a 
diagram showing the field lines often helps us to describe the vector field which 


they illustrate. However the field lines tell us only the direction of the vector field 
and not the magnitude, 


Example 3 
Describe briefly the vector field corresponding to the field lines in Figure 15. 


Solution 

The vector at any point P is along the tangent to the circle through that point (see 
Figure 16), and this direction is perpendicular to the vector OP. Since an 
anticlockwise sense around the circles is indicated by the arrowheads, the direction 
of F is that of —sin@i + cos0j, If P is the point (x,y) then sin® = y/,/x* + y? and 
cos = x/,/x? + y?, so the direction of F is that of (—yi + xj). The vector field 
function has the form F = f(x, y)(—yi+ xj) for some function f(x, y) taking 
positive values only. 


Exercise 6 
Describe briefly the vector field corresponding to the field lines in Figure 17. 
[Solution on p. 44} 


Summary of Section 1 


1. If ¢ is a function which associates a unique scalar to each point in a given 
region of space, then ¢ is called a scalar field function or a scalar field. 


2. A scalar field in two dimensions is defined over a region of the plane and can 
be represented by a function of two variables, (x,y). Such a scalar field can 
be described pictorially by contour curves; these are curves drawn on the x, y- 
plane with equations given by (x,y) = constant. 


3. Avscalar field in three dimensions can be represented by a function of three 
variables, (x,y,z), and can be described pictorially by contour surfaces. These 
are surfaces with equations given by (x, y,z) = constant. 


4. IfF is a function which associates a unique vector with every point in a given 
region of space, then F is called a vector field function or simply a vector field. 


= 


Figure 14 


Figure 15 


Figure 16 


Figure 17 
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5. A vector field can be described pictorially by field lines. A field line is a curve 
whose tangent at each point is in the direction of the vector at the point. 


2 Differentiating scalar fields 


In the study of functions of more than one variable in the last unit, certain 
properties of the function were best understood in terms of the partial derivatives 
of the function. For instance, at a stationary point of a surface the first-order 
partial derivatives are zero. Whether the stationary point is a local maximum, 
local minimum or saddle point depends on some combination of the second-order 
partial derivatives. 


Certain characteristics of scalar and vector fields can be described in terms of 

the spatial rate of change of the fields, and we shall see that these involve a 
combination of the first-order partial derivatives of the field functions. The 
properties that we are going to introduce in this and the next section provide us 
with relationships between a scalar field and a vector field, and between different 
vector fields. These are fundamental properties of fields that will be used in the last 
section and in Unit 30 to characterize and describe the influences of various 
physical quantities in the real world. 


2.1 The gradient of a scalar field 


A scalar field is described pictorially by contours, i.e. contour curves for a two- 
dimensional field and contour surfaces for a three-dimensional field. The 
closeness of the contours tells us quite a lot about how the field is changing in a 
particular direction. For instance, land elevations are represented by contour 
curves drawn on a map. If the contour curves are close together then there is a 
steep hill in that region and if the contours are well separated then the land is 
reasonably flat. 


The change in a scalar field in different directions can also be described by looking 
at the contours. Suppose we consider the elevation of the land above sea level, 
represented by the function (x, y), as an example of a scalar field. Figure 1 shows 
a hill and its representation on a map by contour curves. 


A 
= 
<2” 


x x 


Figure 1 


If we travel round a contour curve from A to B say, then the value of the function 
(x,y) does not change because a contour curve is defined by #(x, y) = constant. 
Now suppose that we put a ball at the point A on the hill. Experience tells us that 
the ball will begin to roll down the hill along the line of steepest slope, in other 
words the direction in which the scalar field is changing most rapidly, We would 
probably expect that this direction is perpendicular to the tangent to the curve AB 
at the point A. The acceleration of the ball will depend on how steep the slope of 
the hill is near the ball’s position. At each point of the hill the magnitude and 
direction of the initial acceleration of a ball defines a vector field. We would expect 
this vector field to depend on how the scalar field of land elevations changes from 
point to point. It is this change in a scalar field that is described by the gradient. 


By ‘spatial’ we mean with 
respect to space coordinates; 
that is, not with respect to 
time. 
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In this subsection we are going to make this idea more precise. We begin by 
adopting the following definition of the gradient: 


At each point in the domain of a two-dimensional scalar field 
we can evaluate the vector 

ap. oO. 

Fra + ay i. 
giving a vector field called the gradient of the scalar field ¢ and 
denoted by grad ¢. 


We will see that the gradient of a scalar field tells us how the scalar field is 
changing. But before we do that, the following examples and exercises show that 
finding the value of grad @ at a point is relatively straightforward. 


Example 1 
Evaluate grad ¢ at the point (0,1) when ¢ = log,(x + 2y). 


Solution 
The vector field grad @ is formed by evaluating the partial derivatives of of first 
order: 
$ = log.(x + 2y); 
a1 a2 
éx  x+2y es @y x+2y" 
1 2 

So grad p = (5) + (; 5 a 


Putting x = 0 and y = 1 we have, at the point (0.1), grad @ = Hi +j. 


Exercise 1 
Evaluate grad @ at the point (—1,2) when ¢ = x?y. 
[Solution on p. 44} 


Example 2 

Evaluate grad @ at the point (1,2) when ¢ = x* + y?. Show that the direction of . 

grad ¢@ at (1,2) is normal to the contour curve through the point (1,2). eae means ‘at right angles 
Solution 


The vector grad @ at (1,2) is found by evaluating the partial derivatives of @ of 
first order at the point (1,2): 


g=x? +57; 
o aa and 2 at (1,2) 
Z- y and ea sat (1,2). a 
So at (1,2) has 
grad = 301+ Oj = + 


Let P be the point (1,2), The contour curve through P (see Figure 2) is a circle of 

radius \/5. The normal to the circle through P is in the direction of the vector 

OP =i+ 2j = 4(2i + 4j), so at (1,2) grad @ is in the same direction as this normal 

vector. Figure 2 
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In fact we will now show that the direction of the gradient of any scalar field @ at dr 
a point is along the normal to the contour through that point. 


Suppose that (x, y) is a scalar field in two dimensions and consider one contour P(x 0.0) 
defined by o (x,y) = C where C is some constant, and suppose that P is a point 
on the curve with coordinates (xo, yo) (see Figure 3). 


contour 


$=C 


In Unit 14 we introduced the idea of a curve in the plane as the representation of 
some vector function of a scalar t; so the position of a point on the contour ¢ = C 
can be written as 


Figure 3 


r=x(t)i+ y(t)j 


where 1 is a scalar, and as ¢ varies the point moves along the contour. A tangent 
vector to the contour curve is then 


With this representation the value of the scalar field (x,y) as we travel along the 
contour curve is also a function of 1, Using the chain rule of Unit 25 we can write 


dp a Og dx ee dy 
dt dx dt tay dt 
which is the dot product of two vectors, i.e. 


a dx, dy) _ de 
= (eet) (Ge Si) 0 


where 


= i + i = grad. 


d di 
Since @ has constant value along the contour, =e = 0; hence Bs =0, and the 


two vectors g and are perpendicular. Thus the direction of grad ¢ is 
perpendicular to the tangent vector, i.e. grad is along the normal to the contour, 


Exercise 2 
Find a vector whose direction is normal to the curve x* — 2xy + y? = 9 at the point (0,3). 


[Solution on p. 44) 


For a three-dimensional scalar field represented by the field function (x,y,z), the 
gradient of ¢ is defined by 


6 
grad = ia Fj +, 


Example 3 

Consider the scalar field (x,y,z) = @ayseye 

where C is constant. Show that the gravitational field F = — Fe, defined in 
Example 2 of Section 1 can be expressed in terms of this scalar field o by the 
relation F = grad ¢. 


Solution 
Cc 0p Cx 
For aa Pape SMES “(t+ y+ zt?” 
a 
a cy ea tes Cz 


éy (+? +22? éz “@+y? +27p?* 
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so that 


Cc hee 
grad d = aoe (xi+ yj + zk). 


+ y¥+2 
This is the vector field function describing the gravitational field at the point 
(9,2). 


This example illustrates an important application of the gradient. Scalar and 
vector field functions that we use to model physical quantities in the real world are 
often related in this way; i.e. some vector field function is just a scalar multiple of 
the gradient of some scalar field function. 


For instance, the gradient of the temperature field on the unevenly heated plate in 
Section 1 is a vector field which determines the flow of heat across the plate. In 
Figure 4 the dotted lines represent the contour curves and the arrows represent 
the gradient vector at various points on the plate. 


Exercise 3 
Evaluate grad ¢ at the point (—1,1,0) for @ = xy + z?x. 
[Solution on p. 44) 


At a point P on a contour surface describing a three-dimensional scalar field, the 
gradient is perpendicular to the tangent plane at P. 


Example 4 


Find the direction of a normal vector to the surface n,x + my + nyz = C at the 
point (1,0,2). 


Solution 

The surface n,x + my + n3z = C is a contour surface of the scalar field 

(x,y, 2) = myx + may + n3z and the gradient of ¢ evaluated at the point (1,0,2) is 
normal to the contour surface. We have 


ti) ti Ci 
pan, ean and Bon at (1,0,2), 


so that grad @ = ni + n2j + nyk. This vector is the same for every point on the 
surface. This result should not surprise you because the equation 
n,xX + my + n3z = C is the equation of a plane, as we saw in Unit 14. 


The gradient of a scalar field ¢ is defined in terms of the partial derivatives of @ of 
first order and its direction is along the normal to the contour curves for a two- 
dimensional scalar field or along the normal to the tangent plane for a three- 
dimensional scalar field. So in some sense it is telling us about the change in # as 
we go along the direction of the normal. 


Consider again a two-dimensional scalar field @; we now investigate the rate of 
change of ¢ in the direction of an arbitrary unit vector e. 


Figure 5 shows several contour curves for this scalar field and at a point P a unit 
vector e. Suppose that the curve C is the path of a particle moving through the 
domain of the scalar field with unit speed, and that the direction of the tangent to 
the curve at the point P is given by the unit vector e, Then the velocity of the 
particle as it passes through P is equal to e. Using the idea of a curve as the 
representation of a vector function of time ¢ as introduced in Unit 14, the position 
of the moving particle on the curve can be written as 


r=x()it+ y(tj 


dr 
and then e = a 


heated corner 
Figure 4 


Figure 5 


= constant 
(contours) 
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Using the chain rule of Unit 25 we have the rate of change of with respect to t, 
at the point occupied by the moving particle, as 
dg _dbdx | ap dy 
dt dxdt dydt 


= (grad g) 4. 


Now grad @ is along the normal to the contour at P, and a is the velocity of the’ 
moving particle which as we have seen is equal to e. 


Suppose that @ is the angle between grad ¢ and e (see Figure 6). Then the aboye 
equation becomes 


“ = (grad ¢)-e = |grad ¢| |e|cos @ = |grad #|cos@ 


since e is a unit vector. 


Considering how this expression varies with @ we see that the rate of change of 
is zero when @ = n/2, ie. along the contour. This is what we would expect because 
the contour is defined by ¢ = constant. The maximum rate of change in # occurs 
when @ = 0, i.e. along the direction of the vector grad #, and this maximum rate 
of change of is |grad |. This means that the gradient of a scalar field is the 
vector field which represents the magnitude and direction of the maximum (spatial) 
rate of change of the scalar field, This maximum rate of change is along the normal 
to (i) the contour curve for a two-dimensional scalar field and (ii) the contour 
surface for a three-dimensional scalar field, 


Exercise 4 
For # = x*y?z? find 


(i) the maximum rate of change of @ at the point (—1,1,1), and the direction in which 
this maximum occurs; 


(ii) the rate of change of @ at the point (2,1, —1) in the direction of the unit yector 
e=i+¢k. 


[Solution on p. 45} 


2.2 The vector operator V 


We now introduce a special notation for the gradient of a scalar field. Instead of 
grad we often write Vp where the symbol V is pronounced ‘del’ or ‘nabla’ and 
é 

oz" 

When written to the left of the symbol representing the scalar field function 

(x, ),2) it means ‘find the first-order partial derivatives of @ and form the vector 


0b | 0h | dp’ 
i Bure tke. 


Bel a, 
Mee at 


Vo 


The outcome of this operation is the gradient of the scalar field function. Usually 
we call V¢ not ‘del ¢’ but ‘grad ¢’. V is not a vector by itself but is called a 
differential vector operator implying that it needs something to ‘operate on and 
differentiate’ in a special way. 


Exercise 5 
We know that if f and g are two functions of a real variable ¢ then 
d df dg 
nae tla 
What would be the analogous result for V(dy) where # and W are two scalar fields? 


[Solution on p. 45] 


grad 


Figure 6 


@ = constant 
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Summary of Section 2 


1, The vector field ey € zi + oy formed from a scalar field ¢ is called the 
gradient of ¢ and this is written as grad ¢ or Vo. 
2. Properties of grad @: 


(i) the direction of grad ¢ at a point P in the domain of the scalar field , is 
either along the normal to the contour curve drawn through P for a two- 
dimensional scalar field or along the normal to the tangent plane of the 
contour surface drawn through P for a three-dimensional scalar field 


(ii) the magnitude of grad @ at a point P is the maximum rate of change of 
at the point P. 


3. The rate of change of a scalar field @ at a point P in the direction of a unit 
vector e is given by (grad ¢)-e. 


3 The scalar line integral 


3.1 What is the scalar line integral? 
You are familiar with the idea of integrating a continuous function over an 
interval [a,b] of the real line, i.e. you can evaluate integrals of the form Tf (x)dx 


which are called definite integrals. The definite integral is defined as the limit of a 
sum in the following way: 


We divide the interval [a,b] into N subintervals and in each subinterval we select 
a point: 
x, in the first subinterval, 


x, in the second subinterval, 


x; in the ith subinterval 
and so on. 
Then we form the sum 

S(&1)O%1 + f2)0%2 +o + [(R IX) Ho + SRW )OXY 
where dx; is the width of the ith subinterval, 


The limit of this sum as we increase the value of N indefinitely is defined to be the 
definite integral of f from a to b and is written 


» N 
i P(x)dx =lim Y f(x))5x;. 
A im 


This section is about generalizing the definite integral to give a rather special 
integral of a vector field function that occurs in applied mathematics; it is called 
the scalar line integral. You will see one application of it in Section 4 of this unit. 


Consider a vector field F and let A and B be two points in the domain of F joined 
by a curve C (see Figure 1), Divide the curve into n segments and suppose P, and 
P,,, are two points on the curve with position vectors r; and r; + dr;, where PP, ; 
is one segment. At each point on the curve C the vector field has a particular 
magnitude and direction; at PF, let the vector field function be F;. The scalar line 


integral of F along C is defined to be lim ¥> (F;-dr;), and denoted by 
1 


f F+dr, or sometimes by F+dr. 
© AB 


The definite integral is defined 
in M101, Block II]; MS283, 
Block 111; TM281, Block IV. See 
also the next unit in this 

course. 


a} 


Figure | 


If F = f(x)i and C is the 
interval [a,b] of the x-axis 
then 


[rede= [soon 
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The curve C is called the path of the integral. The path of a scalar line integral is 
said to be an open path if the two end points do not coincide, as in Figure 2, ora 
closed path if the two end points do coincide, as in Figure 3. 


We can integrate along the closed path in two directions, ie, starting at A and 
ending at B or starting at B and ending at 4. We usually adopt the convention of 
going anticlockwise round closed curves and to remind us of this we write the 


scalar line integral along a closed curve C as f F-dr. 
¢ 


The scalar line integral has many applications. For instance, if a vector field and a 
scalar field @ are such that F = V¢, then given F we can find @ by evaluating the 


scalar line integral i F -dr; we shall be more precise about this later in this 


c 
section. In Section 4 we shall describe one of the properties of a vector field, called 
the curl, in terms of a scalar line integral, 


Before we see how to evaluate scalar line integrals we introduce a method of 
representing every point on a path in terms of a single scalar parameter. 


3.2. Parametrization of curves 


In two dimensions a curve can often be described by an implicit equation 

S(x,y) = 0. For example, f(x,y) = y — 2x — 3 = 0 is the equation of a straight line 
or f(x,y) = x? + y? — 9 = 0 is the equation of a circle of radius 3 and centre the 
origin. In three dimensions this method is not suitable. An implicit equation 
J(%,y,2) = 0 describes a surface, not a curve. 


We shall adopt a method to represent the paths of line integrals which applies to 
curves in two and three dimensions; the method is based on an idea introduced in 
Unit 14. 


In Unit 14 a vector function f of one real variable (t, say) was described pictorially 
by a curve in space so that any point on the curve had a position vector 


r=xi+ yj + zk = f(t), 
or in terms of the components, 

A(t) = f(Oi + Ali + (Dk, 
giving the coordinate equations 

x=fil0, y=hlt, 2 = falt). 


These coordinate equations are called the parametric equations of the curve, and ¢ 
is called the parameter; put in a value of t and we have the position vector of a 
point on the curve. 


A parameter familiar to you from mechanics is time. As a particle moves in space 
subject to certain forces, its position vector is a function of time and the Cartesian 
coordinates of the particle's position at any time serve as the parametric equations 
of the path along which the particle moves. For instance, in a problem of a 
projectile moving in a vertical plane, the position vector of the particle at time ¢ 
after projection is 
1 = (uot + Xo)i + (—4gt? + vot + zo)k 
where (xpi + 2k) and (voi + vk) are the initial position and initial velocity of the 
particle respectively. The path described by the particle has parametric equations 
X=Uott+Xo,, y=0, z=4dgt? + uott+2 


and these equations represent a parabola in the x, z-plane. However, the parameter 
that we use in general need not have any physical significance. 


Example 1 
Consider the parametric equations 


x=2 y=4+1,  z2=2  (tany real number). 


A 
Figure 2; an open path 


Figure 3; a closed path 


18 MST204 26.3 


Find the coordinates of the points for t= —1, —4, 0, 4, 1, and sketch the curve 
that these parametric equations represent. 

Solution 

The coordinates of the points are found by substituting each value of t in turn into 
the equations giving the points (—2,5,2), (—1,2,2), (0, 1,2), (1,2,2), (2,5,2). Every 
point on the curve has z = 2 so that the curve lies completely in the plane z = 2, 
as shown in Figure 4. 


t=-t t=-1 


plane z= 2 


Figure 4 


The parametric equations represent the parabola y = x? + 1 drawn in the plane 
z=2. 


Exercise 1 
Consider the parametric equations 
x=acost, y=asint, z=0, (0<t<2n), 


Find the coordinates of the points for which t = 0, t = 1/4, x/2, 3/4, m, 3n/2, and sketch 
the curve that these parametric equations represent. 


(Solution on p. 45] 


The parametrization of a curve is not unique. For instance we can show that the 
points on a circle of radius a and centre the origin in the x,y-plane can be defined 
by 


(i) x = acost, y =asint (0<t<2z) 
a(l — t?) 


or iim oe 
Ter. 75 aes 


(t any real number). 

In the above examples the parametrization of the curve has been given; but in 
two dimensions a path may be given as an implicit equation and finding a 
parametrization is not always easy. One method that will sometimes work is to 
let x = t (or y = ¢) and then solve the implicit equation f(t, y) = 0 (or f(x, t) = 0) 
for y (or x) as a function of tr, 


For example, the general equation of a straight line is y — ax — b = 0 and if we let 
x=tthen y = at + b. These are the parametric equations of the straight line. 


However, this method does not always give a parametrization of the entire curve. 
For instance, for the circle x? + y* — a? =0, if we let x = t, (for —a <t <a) then 
y’ =a* — #. Solving for y gives y= +,/a? — ¢?. The solution y = +./a? — (7 
gives a parametrization for the upper semicircle (y > 0) and y = — ,/a? — ¢? gives 
a parametrization for the lower semicircle (y <0), but the method does not give a 
parametrization for the whole circle. The parametrization using the trigonometric 
functions, as in Exercise 1, does parametrize the whole circle. 
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Example 2 

The general equation of an ellipse with centre the origin in the x,y-plane is 
2 

_ + i — 1=0. Find a parametrization for this curve. 


Solution 

We could try letting x = at (for —1 <t <1) so that y? = b*(1 — 1°), but then we 
would have to treat the two halves of the ellipse for y > 0 and y < 0 separately. 
For y > 0 we have y = b,/1 — ¢? and for y > 0 we have y = —b,/1 — 17. 

A parametrization which covers the whole curve is to let x = acost (for 

0 <1 < 2n), and then y = bsint. (There are other parametrizations as well.) 


What we are going to do is to use the parametrization of a curve to help us in 
the evaluation of a scalar line integral. Any parametrization will do, but for some 
the resulting integration is much easier than for others, 


3.3 Evaluating the scalar line integral (Tape Subsection) 


On the audio-tape that follows, we shall see how to evaluate the scalar line integral, 


f F + dr, for different vector fields F and different paths C. But first we shall show 
Cc 


that by representing the path in terms of a parameter t, the scalar line integral can 
be written as a familiar looking integral in terms of t. 


We defined the scalar line integral as the limit of a sum of the terms F,- 6dr, as the 
number of line segments gets larger and larger. The position vector r, of the point 
P, on the path C (see Figure 5) can be given in terms of a parameter t by 


Hy = H(t)) = x(t + ye) J + z(t) k, 
where ¢; is the value of t corresponding to the point PF. 


As we move along the curve, t changes in value and the change in r, for a small 
change in ¢ is just 


br = sri + Bony + Stem 


and we can write 
sf ! dr, 
YFron= ¥ (r #) 5. 
fet fen 
In the limit as n increases indefinitely we can write 


tm (SE) = f° (ret) a 


where t, and f, are the values of at the points A and B respectively, 
The integral on the right hand side can be converted to a familiar looking integral; 


if we let ret = f(t) then the scalar line integral becomes 


li S(t)dt. 


The audio-tape develops an algorithm for evaluating scalar line integrals; it shows 
how to make the transformation from 


[ra to iy (r-# dt. 


Start the tape when you are ready 


Figure 5 


dx, dx 
By7 we mean rr evaluated 
at ,, that is, evaluated at “4 
point P,, and similarly for z 
dz, 


id 
at 
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o A Tiwo- dimensional line integral 


Evaluate : E.dr 
Cc, 


where F = (2x+y)i-xj 


and ¢, is the straight line befiween (2,0) and (0,2) 


y= ox +b 


a= gradient = 


b= y intercept = 


A 
@,0) * 


equation of line is 
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gradient a=-l 
y intercept b= 2 


equation y=a-x 


oO Parametrize 


y=2-x 

x= 

y= 

G; goes from (2,0) To (0,2) 


C goes From | to 
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ts) Substitution 


at 
Finally t= 2 and t,= O 


, LJ 
Hence [e-s)s Me [et Saal a 


Solution at the end of the Frames 


MST204 26.3 Tape 23 


Ce) Evaluation 


fe +2t) dt 


look af path 


parametrize of integration 


write all the 


subetitote variables x,y (etc) 
in terms of parameters 


this looks more like 


‘ evaluate what I’m used to 


(3) Another 2-d problem 


Evaluate i E.dr 
¢, 


2 


where F =(2xt+y)i - xj 


ond C, is the quarter circle, centre at 
the origin and radius 2, beliseen 


(2,0) and €,2) 
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Equation xt + y* =4 


y= 


Cz starts af (2,0) and 
ends, at (0,2) 
A Sot goes From 


© Substitution 


t, 4 
ie (Ge st) at 


Qx+y)i —xj 
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Solution at the end of the Frames 


(2) Comparison of answers 


BA 
(0,2) L® 


© Another vector Field 
evaluate j, Fide 
a 


where F= x? yj 


a) for C =C, di) For C=C, 
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2 2 
r)at =| (t+ t-2) de=-F 
2 


— 


y= 2sint 
EB = 4 costiti+ 2cin tj 


c= -2sintit2 costj 
tb 1/2 

dr oa 2 . . 2 
F.2 )dt =|(-8 tsint +4 cost sint) dt =-2 
(F ¢) Js cos in OS int) 5 


u ~ 
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$ E.dr=? 
Cc 


where F = xye + 2) 


and C is the closed curve 


Hh { Erdy: 
APB 


Parametrize Parametrize 


(0.2) (2) 


(2,0) 
substitute substitute 
t. 
f"€. 2)at 
t, dt 


evaluate 


ty 


evaluate 


J F.dr + 
APB 


BQA 
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Answers 


i) f eatels 
APB 


Parametrize 
fe x 2 cos t 
e 
2 sin t 
3 O and t, = "% 
substitute 


= 4 sin Ecostit2j 


—2sinti +2 cost) 
evaluate = 


TW 5 
f (-8 sin? t cos t +4cos Eat = 3 


Porametrize 
Yte 
Q 
A 
substitute ~ 


evaluate 
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What happens for Fextit+y) & 


for every 


CONSERVATIVE FIELD 
closed path 


For some 
closed path NON-CONSERVATIVE FIELD 


Summary : How To evaluate line integrals 


step | 2 parametrize 


step 2: substitute 


step 3: evaluate 


This applies to three—dimensional line integrals 


as well as two-dimensional line integrals 
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Solutions to Exercises in the tape 
Frame 5 
F=(2+i-tj 
r=ti+(2—8j 
dr a 
a 
Feta et oti 2+2) 
=2 and =0 


0 
f (2 + 2n)dt 
2 


Frame 11 


fics —8sintcost)dt [-« —4sin? iT" 

= (-2n-4)- (0) 
= —4-2n. 

Frame 15 


The path APB is the same as the path C; and Frame 13a gives the solution to the line 
integral of F = x*i + yj along C,. We have 


foFede= fo Pede= 3. 
cs Ary 
The path BQA is the same as the path C, but here we go along the path from B to A. In 


Frame 13a we integrated from A to B and found that | F+dr= Fedr = —4. 
ia) ‘AQ8 


The integral from B to A is just minus this answer because we just reverse the limits. 


Thus | F-dr=4. 
BOA 


34 Post-tape exercises 
Exercise 2 


Evaluate the scalar line integral f u-dr for each of the vector field functions u given below 
Jc 
and where C is a circle defined by the parametric equations 
xX =acost, y =asint, z=0 (O<t<2n). 


(i) w= (xi + yi); 


aes 
x+y 
(ii) w= —yi + x. 
{Solution on p. 45] 


The solution to this exercise illustrates one important point that was introduced 
on the tape: that the value of a line integral round a closed curve may not be zero. 


This exercise and the problems in the tape section are examples of line integrals 
along paths drawn in the x,y-plane — these we called two-dimensional line 
integrals. We shall now look at some paths in three dimensions. The same three- 
step approach that we introduced on the tape applies here; the steps are 

(i) parametrize, (ii) substitute and (iii) evaluate. 


Example 3 
Evaluate the scalar line integral of F = xyi + yz*j + y?zk along the path C 
between the points (0,0,0) and (1,1,1) defined by the parametrization 


x=t y=r, z=t (O<r<1). 
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Solution 


The parametrization of the path is given in the question; this is usually the case 
for curves in three dimensions —- they are defined in terms of a parameter. The 
limits of the integration are t, = 0 and f = 1. 


Now we must replace x, y and z throughout in terms of t. We have 
F=fit+j+ek 
and or=xityj+zk=ti+ Pj +cek 


so that 
dr ~ 
a lt 2j+k 


di 
and Peo=e +28 +8 = 2 4 325, 


On substitution the line integral f F+dr becomes 
c 


r 
f (8 + 30)dt. 
0 
This is a familiar looking integral and evaluating it in the usual way we get 
1 
it F-dr= f (C+ 3¢5)dt = 3. 
¢ 0 


Exercise 3 

Evaluate the scalar line integral of 
F = (2x)i + (xz — 2)j + xyk 

along the path C, between the points (0,0,0) and (1, 1,1) defined by the parametrization 
x= y=, z=P O<t<t) 

[Solution on p. 45) 


Exercise 4 


In Exercise 3 the path of the integration starts at (0,0,0) and ends at (1,1, 1). Repeat the 
exercise for the path starting at (1,1,1) and ending at (0,0,0) and following the path of 
Exercise 3 in the opposite direction. 


[Solution on p. 45) 
The solution to Exercise 4 illustrates an important property of the scalar line 


integral. If we reverse the end points of the path then the value of the line integral 
changes sign, ie. for the path C in Figure 6 


Exercise 5 
(i) For the vector field defined in Exercise 3, evaluate f F-dr for the path C, between 


ce 
the points (0,0,0) and (1,1,1) defined by the parametric equations 
x=h yeh z=P (<r<1) 
(ii) Use the results of Exercise 3 and part (i) to find the value ot F-dr where C is the 


ic 
closed curve starting at (0,0,0) going along C, to (1,1,1) and returning to (0,0,0) 
along the path C;. 


[Solution on p. 45] 


Example 4 

Show that the scalar line integral f F -dr for the vector field F = x7i + yj is 
AB 

independent of the path taken between A and B. 
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Solution 


Consider two points A and B with coordinates (x,,y,,z,) and (x2,)2,23) 
respectively, and any curve C joining A to B (see Figure 7). 


Let a parametrization of the curve C be 
x=f(), y=alt), z=h(e), 
with ¢ =f, and t = t, at A and B respectively. 
Then F = (f())Pi + g()j, 5 
. dr d, de, ah Figure 7 
dt dt dt” dt 


2 de 
s0 that f F-dr = f Fedr=[ (Fe) at 
AB c in \ at 
ef of | dg 
= (law. 8 
-[ (7 at Fae] 


Slee: 
-[F+8] 


A(X Iz) 


_ (Ue) , (ee)? (MeoyP (gt)? 
~ (A, eco) _ (er, tet) 


Now /(t;) and f(t2) are the x-coordinates of the points A and B; also g(t,) and 
g(t2) are the y-coordinates of the points A and B. The value of the line integral can 
thus be written in terms of the coordinates of A and B as follows: 


3 2 3 2 
%2 V2 bail VN 
Pn (ec ese REST PG 
bts’ ti 344) (f+4 


This result depends only on the positions of the end points A and B and not on 
the path joining A and B. 


The vector field in Exercise 5 is an example of a non-conservative vector field 
because in part (ii) we found that the value of its line integral round a closed curve 
is non-zero, Compare this with the vector field defined in Example 4. Consider any e 

two paths labelled APB and AQB joining A to B (see Figure 8). The value of the Figure 8 
line integral along each of these two paths is the same: as we found in Example 4 

it is 


x2 yi) [xt ot 
a cae ia Eh 


Hence the value of the line integral of x?i + yj round the closed curve APBQA is 
zero because 


f Fedr= F-dr+ [ Fedr= f Fedr— F-dr =0, 
APBQA APB BOA APB AQB 


and this will be true for every closed path. The vector field F = xi + yj is an 
example of a conservative vector field. 


Exercise 6 
(i) Show that the vector field grad @ where ¢ = 


2 
“ 5 is the same as the vector field 


3 


given in Example 4. 


(ii) Show that for the vector field grad @ in part (i) and the points A = (x;,),,2,) and 
B= (X2,)2,22), 


J wrad 6: de = d(e2,92520) — blr. 20) 
he 


[Solution on p. 46] 
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Exercise 7 


Use the same method as in Example 4 to show that the value of the line integral of the 
vector field grad @ for any scalar field is given by @(B) — (A), where A and B are the two 
end points of any path in the domain of @. 


[Solution on p. 46] 


The result of Exercise 7 shows that the vector field grad @ obtained from a scalar 
field @ is a conservative vector field. We shall continue the discussion of 
conservative vector fields in Subsection 4.3 of this unit. 


Summary of Section 3 


1. The scalar line integral is of the form f F + dr where F is a yector field function 
and C is a curve along which the integration is carried out. If C is a closed 
curve then we write $e +dr. 

2. A curve can be expressed in terms of a parameter rf by the parametric 
equations 

x=x(0, y=yt) z=2(t) (a<t<b) 

3. To evaluate a scalar line integral we express F and r in terms of a parameter ¢ 

and then the line integral rede can be written in terms of an ordinary 


integral in ¢, as summarized in the following procedure: 


Procedure 3.4; Evaluating line integrals 


Step 1. Parametrize: find a parametrization of the curve so that 
every point on the curve is given in terms of a parameter 
t 


Step 2. Substitute: write each variable x, y, z, the vector field F 
dr 
dt 
an integral in terms of t. 


and — in terms of ¢ and write the line integral freae as 


Step 3, Evaluate: evaluate the ¢ integral. 


4. The line integral of the gradient of a scalar field @ is independent of the path. 
We have f Vo-dr = $(B) ~ $(A), 
AB 


5. A vector field F is conservative if for every closed path C drawn in the domain 
of F we have 


Fede =o. 
é 


Otherwise F is a non-conservative field. 


4 The curl of a vector field 


4.1 Line integrals and curl (Television Subsection) 


The television programme investigates line integrals of velocity fields representing 
three fluid flows and uses the line integral of a vector field u to define a new vector 
field called curl u. 


During the television programme the results of the following exercises are used; so 
work through the exercises before watching the programme. 
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Exercise 1 » 
A circular disc rotates about its centre O with constant angular velocity «; that is to say, 

the angle turned through in time ¢ is wt (radians). Any point P on the disc has velocity of 
magnitude ma and in a direction perpendicular to the radius vector OP where a is the 

magnitude of OP (see Figure 1). P 


wa 


(i) Show that the velocity vector field u is given by u = @(—yi + xj) using Cartesian 
coordinates with origin at 0. 


(ii) Evaluate the line integral u-dr where C is a circle whose centre is the origin and Oo 
€ 7 
whose radius is b. Figure 1 


[Solution on p. 46] 


Exercise 2 
If f(y) is a differentiable function of y, show, by using Taylor series, that 
oy’ oy If you find this exercise 
F ily = * =(y Dicey ao difficult, consult the solution. 
im _—1_1 _“! 2. 
ay0 by dy 


[Solution on p. 46] 


The programme can be summarized as follows. 
Part 1; we begin by introducing two different flows of a liquid: 


(i) the motion of water in a bow! which is rotating with constant angular 
velocity, and 


(ii) the motion of water in a vortex tank, which is rather like what happens to 
bath water near to a plughole. 


These two flows appear similar because the motion is circular; but small floats 
placed in the fluid behave differently in the two cases. 


Part 2: the rotation of a rigid disc is used as a simple model for the motion in the 
bowl and we find the velocity field u representing this flow. We then use a line 
integral to define the z component of a new vector field called curlu. The 
expression defining (curlu), relates the line integral to the rotational properties of 
the vector field u. The value of (curlu), is used to predict whether small floats 
placed in the fluid flows will rotate or not. 


Part 3: the flow in a long rectangular tank is then investigated. At first sight there 
appears to be no rotational motion of the fluid in the tank and we would expect a 
float not to rotate. However the formula for (curl u), predicts otherwise. Have we 
made a mistake? 


Now watch the television programme and find out, 


In te television programme we defined the z component of the curl of a two- 
dimensional vector field in the x,y-plane. We can define the curl for any such 
vector field in the following way: 


TV26 


Suppose that F is a two-dimensional vector field defined over a 
region of the x,y-plane. Then the z component of the vector field 
curlF at a point P in the domain of F, is defined by 


aie 
(curlF), = lim (br a 


where C is a curve drawn in the domain of F surrounding the 
point P and A is the area within C (see Figure 2), “ GQ 
4 & 


The result is the same regardless of the shape of the curve C. 


Figure 2 


During the programme the results of several calculations involving the limit in this 
definition above were assumed. In the following exercises you are asked to show 
that these assumptions are correct. 
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Exercise 3 r 
The velocity field for a disc rotating about the origin with constant angular velocity « is a 
1 D 
given by w= — wyi + wx}. Show that lim a was) = 2 in each of the following cases: 
aol AS 
(i) Cis a square of side b, and centred at the origin (see Figure 3); 
(ii) C is a square of side b, and centred at a point (xo, yo) (see Figure 4), 7) ‘por 
[Solution on p. 46] 
Exercise 4 A B 
‘The velocity vector field for the flow in the vortex tank is given by ———— 
: Figure 3 
u= a” u 
D c 
where e, is a unit vector perpendicular to the radius vector at a point P (see Figure 5) Ms 
distance d from the origin. ' y 
rs 
(i) os 6 
(xo.¥0) 
(ii) Hence show that u=—; 
(+) 

(iii) Evaluate the line integral of u round the curve C = ABCD shown in Figure 6, where A B 

AB and CD are parts of radii and BC and DA are arcs of circles, 0 
(iv) Hence evaluate (curlu), at the point P. Figure 4 
[Solution on p. 47] vp 

un Sey 


We saw in the television programme that the curl of a velocity field u describes the 
rotational properties of the field. We could predict whether a small float placed in UN Pcxy) 
a fluid would rotate or not depending on whether the value of (curlu), was non- " 
zero or zero respectively. 


The mathematical analysis for the third fluid experiment, involving the linear tank, 


introduced a different way of calculating curl. We found that for the velocity field. : 
d ‘ Figure 5 
u =f(y)i, the z component of curlu could be expressed as 2 i.e. in terms of the C 
C m yh 
rate of change of the velocity field. 
B 
D 
/ ZA 
(ra 
0 x 
Figure 6 


Calculating derivatives is certainly easier than evaluating line integrals and taking 
limits. In the following example we obtain a formula for (curlF), for the more 
general two-dimensional vector field F = F,(x,y)i + F(x, y)j, in terms of the 
spatial rates of change of the components of F. 


Example 1 


Suppose that F is a two-dimensional vector field defined over a region of the x, y- 
plane. By integrating round a small rectangle show that 
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é 
IR) 
(curl F). tx 
where F, and F, are the x and y components of the vector field. 


Solution 


To find the z component of curlF we need to evaluate lim 5 F-dr for some 
Aw0 ic 
area A contained within a closed curve C. Consider a two-dimensional vector field 
F defined over a region of the x,y-plane and let F = F,i + F,j. For the curve C, 
consider a small rectangle ABCD with sides h and | and centre at the point P with 
coordinates (xo, Yo) (see Figure 7). To evaluate the line integral of F round the 
rectangle we shall appeal to the definition of a line integral as the sum of terms of 
the form F+ dr where dr is a small displacement along the curve of integration. For — Figure 7 
the rectangle ABCD, since we are assuming that h and y Jare small, Il, suppose that we 


divide the rectangle into four small displacements AB, BC, CD, DA where 


AB = hi, BC = lj, CD = —hi and DA = —lj, 


We now approximate the line integral as the sum of four scalar products in the 
following way: 


$ F+dr = Fy" AB + Fyc*BC + Fep* CD + Fy4*DA 
c 


where F442, Fac, Fep and Fp, are the values of F at the mid-points of the sides of 
the rectangle. 


Now 

F45*AB =Fyg*hi= Flaws - i) 
because F,,*i is the x-component of F evaluated at the mid-point of AB, which 
has coordinates (x07 = i) 


Similarly Fyc* BC = [0+ jo 
_ i 
Fep*CD = =A (sour + i) 


— h 
F,4°DA = =F(v0 - ft 


and so we have 


I i 
Fede = (F(x. =) ;) = Flore + Al 
h h 
+ (1[ + jo] = Fe = fava} 


The area contained within C is the area of the rectangle ABCD and this is hl. Thus 


u if 
Fi{X0,Yo — 5] — Fil Xoo +5 
a4 Tee aes ea) |S el ea 

AJ a ij 

h h 
F{. ae 3] — a Xo — 30] 
h 
Now in the limit as A—+0, both h—+0 and |—+0; and we have 


4 
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I u 
Fi(xou i ;) = Filo ch } 


hilo Ui oy These results follow either by 
expanding the numerator as a 
Taylor series about (xo, yo) 


and and then taking the limit as 
h h hy +0, or from the 
= ae efinition of parti 
; (s+ 5.) B(% 3) aF, differentiation in Unit 25, 
lim ===, 
er h 6x 
1 6: 
Finally we can write lim ab reae = ah + es and so 
An0A dy ax 


for a two-dimensional vector field. 


é. 
Note that if F, = 0 then (curlF), = > which is in agreement with the result 


obtained in the television programme for the velocity field of the flow in the linear 
tank, 


We now have a formula for the z component of curlF in terms of the derivatives of 
the components of F. In the next exercise you can check the results of Exercises 3 
and 4 by using this derivative formula. 


Exercise 5 
Use the derivative formula to find the value of (curl u), for the velocity vector fields 


(i) w= —oyi + wxj; 


a ic 
(i) a= way! + xj). 
In each case compare your answer with the relevant result obtained in Exercise 3 or 4, 


[Solution on p. 47] 


The velocity field given in the first part of this exercise represents the motion of 
the rotating disc. We find that (curl u), = 2w at every point on the disc, and that is 
the result we obtained in the television programme by using line integrals. The 
velocity field given in the second part of this exercise represents the flow in the 
vortex tank, We find that (curl u), = 0 in this case, which again agrees with the 
line integral calculation. 


So the curl of a vector field describes the rotational properties of the vector field. 
Its effect can be illustrated by placing small floats into a moving liquid. Then the 
curl of the velocity field at a particular point is twice the angular velocity of a float 
placed at that point. 


4.2 Curl of a three-dimensional vector field 


The television programme and the definition of curl in the last subsection 
considered only two-dimensional vector fields defined over regions of the x,y-plane, 
We defined the z component of curl for a vector field F = F,i+ Fj by the line 


integral k 
1 y 
lim a F-dr, 
a0 A Je * 
Now the z direction is normal to the area A (see Figure 8); and the unit vector k Figure 8 


points in the z direction so we could have written (curlF), as (curlF)-k. 


We define curl for a three-dimensional vector field in an analogous way. 
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Suppose that F is a three-dimensional vector field defined over a 
region of space. Then there exists a vector field curlF with the 
following property: for any point P and any unit vector n, the 
component of curlF in the direction of n is given by 


(curlF), = (curlF)-n = tim > F-dr 
Gote 


where C is a small curve lying in a plane and enclosing the point eS 
P such that n js normal to this plane containing C and A is the Figure 9 


area of the plane enclosed within C (see Figure 9). ‘Wheniale themnie-weelor 


then (curlF), is just what we 
denoted previously by 

This is rather a formal definition which is a simple extension of the definition of (curlF),. 

curl for a two-dimensional vector field. You are more familiar with vector fields 

being given in terms of a vector function and its domain. The domain of curlF is 

the same as the domain of F because curl is defined for the same set of points as is 

F. To write curlF as a vector function of the form a, i+ a,j + a3k we need the 

components of curl in the directions of i, j and k. Now we can find these 

components by using the definition above: 


a, = (curlF)-i= lim a F-dr, 
4,704; Jey 


a, = (curlF)-j = lim + Frdr 


A042 Sy 
and a; = (curlF)*k = lim ak Fedr 
As As dey 


where C,, C> and C3 are curves drawn parallel to the y,z-plane, x,2-plane and x,y- 
plane respectively, and A,, A) and Aj are the areas contained within C,, C, and 
Cy respectively (see Figure 10). 

Figure 10 = z 3 


So we can write 


1 | 
curt = (im, eas} + (tim > Ft; 
a0 Ade, a0 Azle, 
1 
+ | lim —— Fa] k. 
Ay0 Az Jc, 

For a two-dimensional vector field we found that we could write the z component 
of curl in terms of derivatives. Now we show that we can do the same here, 


Example 2 
By taking a rectangle parallel to the x,y-plane show that 
“ee oF, OF, 
lim —® F-dr =—— -—— 
Ast0 As) cy Tox Oy 


where F, and F; are the i and j components of F respectively. 


Solution 
Consider a vector field F and let 


F=FKi+hj+ hk. 
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For the curve C; take a small rectangle ABCD with sides h and | and centre at the 
point P with coordinates (xo, ¥o,Zo), the rectangle being drawn in the plane 
2 = Zo = constant (see Figure 11). 


The analysis is similar to that used in the solution to Example 1, We write the line 


integral () F+dr as the sum of four scalar products as follows: 
Gs 


F-dr ~Fyg* AB + Fac* BC + Fep*CD + Fp," DA. 


Cy 


For this vector field F we have three components; however the four scalar 
products involve multiplying by either i or j so that the k component of F, ie. F;, 
does not appear in the value of the line integral, For example, 


ae 4 I 
Fg" AB =F ap: (hi) = Filxoue = ea) 


and Pye" BC =F gc (Ij) = {re + Foyezo)h 


! 
(v0 -4 


U 
(xs. * i) 


Thus 4 Fedr~= 


adey ! 


h h 
Fe 3 ete) = F(x = 3.02) 


be h 
and in the limit as A, —+0 (i.e. h——+0 and !—-+0) we have 
coe OF, OF, 
lim — Fedr = <2 - <1 
Pee oe wy ox dy 


and this is the k component of curlF. In a similar way we can obtain derivative 

forms for the i and j components of curlF, For the i component we choose for C, 

a small rectangle in the x = xo plane (see Figure 12). The i component becomes 

e - &. For the j component we choose for C, a small rectangle in the y = yo 
OF, OF, 


plane (see Figure 13). The j component becomes 7 er 


Thus we can write the vector field curlF in a form that is easy to evaluate, namely: 


OF, OF; OF, _ OFy 
dz ox}! 


af ah), 
éx ay] 


Applying this formula to the two-dimensional vector field F = F,(x,y)i + Fa(x, y)j 


OF, 36 
we obtain the same result as in the last subsection; that is (curlF). = a - ze 
but now we can go one stage further and write curlF as the vector 

OF, OF 
IF = (|<? 1) x. 
curl | Be aay k. 


We were not able to do this before because we knew nothing about the x and y 
components of curlF. However we could have made an educated guess at their 
values. We saw in the television programme that for the rotating disc, 

(curlu), = 2m. Now angular velocity is a vector quantity; for the rotation of the 
x, y-plane the angular velocity vector is wk, so we might have suspected that 
curlu is also in the k direction. 


Figure 13 
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Example 3 

If F = (xy + 27)i + x?j + (xz — 2)k 

(i) find cirlF at the point (1, —1,3); 

(ii) find points where curlF = 0. 

Solution 

We use the differential formula to calculate curlF. 
R=axyt+?, Rh=x*7, KRex-2, 


OF, _ OF, | 
so that ie id ye 2 
oF, _ oF, _ 
oie 2x and om 0, 
oF, OF, _ 
Gx 77 and aoe 


giving curlF = (0)i + (2z — z)j + (2x — x)k. 
(i) At the point (1, —1,3) we have curlF = 3j + k. 
(ii) curlF is zero when x = 0 and z = 0 for any value of y. 


Exercise 6 

Find curlF if 

(i) Fer; 

(ii) F = s(r)ts 

(iii) F = (y? + 22)1 + (xy + 62)j + (2? + 2xz + yk; 

(iv) F=grad@ where ¢ is a scalar field, 

Here r= xi+ yj +zk, r= |r| = Jttyee and # is a unit vector in the direction of r, 
[Solution on p. 48} 

The type of vector field function in part (ii) of this exercise can be used to 


represent many physical quantities in nature. For example, you will see in a later 
GMm. 


unit that the gravitational force field can be written in the form F, = — aR F so 
that f(r) = — oe. The result of this exercise shows that curlF, = 0. 


We will see in the next subsection that if the curl of a vector field is zero for all 
points in its domain, then the vector field is conservative, so we will use curl 
mainly as a test for conservative fields. 


The derivative formula for curlF contains terms involving partial derivatives. The 
expression is not particularly easy to remember; however the following exercise 
shows an alternative notation in terms of the differential operator del. 

Exercise 7 

By expanding the cross product of V and F show that curlF = V x F_ 

[Solution on p. 48) 

This is an alternative definition for the curl of a vector field which you may come 


across in other texts. The way we have defined the curl in this section more readily 
shows its important physical significance. 


4.3 Conservative fields 


ot is one that we have considered at various 


times in this unit. In Example 2 of Section 1 we saw that this vector field function 


The vector field F = 
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represents the gravitational force on a body a distance r from the Earth’s centre. 
In Example 3 of Section 2 we showed that it could be derived from the scalar field 
o= aa In Exercise 6 of Section 3 we saw that the line integral of F 
V+ y+ 2? 
between two points A and B is independent of the path taken. And Exercise S(ii) 
of this section shows that curlF = 0. Now we shall put these ideas together by 
talking about a special class of vector fields called conservative fields for which 
these last three properties are true. 


On the audio-tape we defined a conservative vector field F as one for which the line 
integral F-dr =0 for every closed curve C drawn in the domain of F. (The 
(e 


reason for the name has to do with the conservation of energy for bodies acted on 
by forces which are modelled by conservative vector fields.) If a vector field does 
not have this property then it is non-conservative. The velocity field u representing 
the rotating disc is an example of a non-conservative vector field, because when C 
is a circle of radius b and centre the origin, 


$ usdr=2onb? —_ (see Exercise 1), 
¢ 


This definition of a conservative vector field has an important consequence for the 
curl of F, We have seen that 


1 1 
1F = tim (2 Fedr}i+ hi a Fat) 
or tin (7 cy | sin (3, a 


1 
+ lim (if Fa] k. 
Ay 0\ Ax Jey 
Now for a conservative vector field, the line integral round every closed curve is 


zero. This means that each of the components of curlF is zero, and we have the 
important result: 


If F is a conservative vector field then curlF = 0, 


The converse of this is also true: if curlF = 0 and the domain of F is sufficiently We shall not say precisely 

simple then F is conservative; however we will not be able to prove this statement what we mean by ‘sufficiently 

in thi This gives an alternative definition of a conservative vector field pal eee con fal Bs 

in this course. This gi a domain we consider in this 

which you may come across in other texts. unit is sufficiently simple in 
this context. 


If curlF = 0 at every point of the domain of F, then F is said to 
be a conservative vector field. 


Thus if we want to decide whether a vector field F is conservative or not, instead 
of evaluating line integrals round every possible closed curve that could be drawn 
in the domain of the vector field, we evaluate curlF in its derivative form and see if 
its value is zero or not. 


Example 4 
Which of the following vector fields is conservative? 
(i) Fs 2x7 + 2y7j +409 + dk; 
(ii) G=xit xyj. 
Solution 
We have to evaluate curlF and curlG. Using the derivative formula for curl gives: 
(i) curlF = (y? — y?)i + (x? — x?)j + @—O)k 
=0i + Oj + Ok =0. 


Thus curl F = 0 so that F is a conservative vector field; 
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(ii) curl = (0 — 0)i + (0 —0)j + 2xy — 1k 
= (2xy — 1k. 


Thus curlG # 0 at every point, so G is a non-conservative vector field. 


Exercise 8 

(i) The vector field F = 2xi + (xz — 2)j + xyk was shown to be non-conservative in 
Section 3 (Exercise 3) by evaluating the line integral along a closed path. Check this 
result by finding the value of curlF. 


(ii) Is the vector field G =r"? conservative, where r = xi + yj+zk, r= |r| and ® =r/r? 


[Solution on p. 48) 


Exercise 9 
If F is a conservative vector field show that the line integral of F between two points A and 


Bie. f F dr, is independent of the path taken, 
Jan 
[Solution on p. 48} 


Exercise 9 shows that for a conservative vector field the value of the line 
integral f F- dr depends only on the end points A and B of the path and not on 
AB 
the path itself. In particular we can define a scalar field at every point 
P = (x,y,z) in the domain of the conservative field F by letting 
(x,y,2) = i) Fed (eee Figure14): 
op 

(1) 
Then { F-dr= $(B) — (4). 

AB 


We recall from the result of Exercise 7 in Section 3 that the value of the line 
integral of grad along any curve between two points A and B is @(B) — @(A). 
That is, 


J, erad ddr = $B) ~ $14), 
AB 
In general it can be shown that when Equations (1) hold then F = grad @ for the 


scalar field ¢. 


One application of this idea is to the motion of a particle acted upon by a 
position-dependent force which can be modelled by a conservative vector field. 
Such a vector field is called a conservative force field. For instance, consider a 
body moving in a one-dimensional force field F = F(x)i. Any one-dimensional 
force field is conservative and the scalar field which we define from this force field 


is just $(x) = f Fows In Unit 4 we defined the potential energy of a body to be 
0 


U(x) = -— f F(s)ds, so the one-dimensional scalar field that we have introduced 
here is just ine the potential energy of mechanics. 

We can generalize the ideas of one-dimensional mechanics and we use the line 
integral — freae to be the potential energy U of a body moving in a conservative 


force field F. Conversely, given a potential energy function U(x, y,2), the force 
acting on a particle is 


F=-VU. 


In a later unit we will follow this idea up. You will see that for a conservative force 
field we have a conservation of energy equation, 


kinetic energy + potential energy = constant, 


and thus the name ‘conservative’, 


* 
Figure 14 
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Summary of Section 4 


1. For a two-dimensional vector field F defined over a region of the x,y-plane, 
the z component of curlF at a point P is defined by 


(curlF), = fim (Gf Fa 
a0 Ate 
where C is a small curve drawn in the domain of F surrounding the point P 
and A is the area within C. If F = F,i + F,j then an alternative way of 
writing (curlF), is 
OF, OF, 
(curlF), = im aie 
2. For a three-dimensional vector field F we define the component of curlF at a 
point P in the direction of a unit vector n to be 


A40 


1 
(curlF), = (curlF)-n = tim( 4 Fa 
c 
where C is a small planar curve enclosing the point P, A is the area enclosed 
within C and n is normal to the plane of the curve C. 


3. We evaluate curl of a vector field by using the derivative formula: if 
F=Ri+hj+Kk 


a. 7} é. 
then eu = Vx F = (SE), (28 Bi (2 Fy 


dz ox}! * Vax ~ ay 


4. The curl of a vector field F is another vector field defined over the same 
domain as F; curl is a detector for rotational motion, For instance, the curl of 
a fluid velocity field describes the spin of a float placed in the fluid. 


5. For a conservative vector field: 


(i) curlF=0 (and conversely, if the domain of F is sufficiently simple); 
(ii) if F+dr = o(B) — (A) where F = grad ¢. 
AB 


The curl of a vector field F can be used as a test to see if F is conservative, 


hE 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. The temperature T at a point r = xi + yj on the metal 
sheet is defined by Tx log, r. But T=0 when r = I, and 
T = 100 when r = 4. Putting T= a(log,r) — b we get b =0 
and a = 100/log,4. Thus the scalar field function is 

_ 100log.r 

~~ Tog.4 * 
Since r = ,/x? + y*, the scalar field is defined over the region 
of the plane given by the set of points (x,y) such that 


is/+y <4 


2. The contour curves are given by 


xy = constant = c 


Fi € 
ie y= 
x 


ay 


The figure shows the contour curves for ¢ = $,¢ = —4 and 
c=. 


3. (ia) There is a hill on the left, near ‘contour 100’. 


(b) There may be a lake where the contours are marked with 
negative numbers indicating land levels beneath sea level, 


(ii)The land should be fairly level on the right; the contour 
curves join places at levels 0, | and 2 over quite a large area. 


(iii) If we start off at point A and walk in the direction of the 
arrow the land begins to climb very quickly; the contours are 
close together so the heights change very quickly indicating a 
hard climb! 


4. The vector quantity is the velocity of the water at points 
on the plate. The vector field function u (ie. the velocity 
field) is defined for points P on the plate, The region over 
which the vector field is defined is the set of points (x, y) 
inside the circle x? + y? = 100 and excluding the origin, The 


yi 


velocity vector field u is directed along OP and has 
magnitude C/|OP| for some constant C. So 


COP c 


(xi + yj). 

5. The Earth’s gravitational force field is directed towards 
the centre of the earth, The left-hand figure shows arrows’ 
representing the gravitational force at various points. The 
field lines are shown in the right-hand figure. The 
arrowheads indicate the direction of the force field. 


6. If you look back at Exercise 4 the field lines in this 
exercise would describe the velocity vector field of that 
exercise. So the direction of the vector at each point is in the 
direction of xi + yj. The vector field function has the form 
F = f(x,y)(xi + yj) where /(x, y) is some function of x and y 
taking positive values only. 


Solutions to the exercises in Section 2 


1 
Ob oxy and Bed at(—1,2); 
Ox 
oF a and oe at (—1,2). 
Vv oy 
So 
2%, ob, 
gadpa s+ a i= 4+) at (—1,2) 


2, The curve x? — 2xy + y? = 9 is a contour of the scalar 
field @ = x? — 2xy + y?, The vector grad @ evaluated at (0,3) 
is normal to this scalar field; 


grad @ = (2x — 2y)i + (—2x + 2y)j 


and at (0,3), grad = —6i + 6j, and this is the vector 
required, 


3. p=xty+ 27x; 


oe 
ae 2xy+27 and 


as at (—1,1,0); 
and at (—1, 1,0); 
and at (—1, 1,0), 
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So grad @ = —2i+j + 0k at (—1,1,0). 
4, (i) grado = 2xy*z*i + 2x*yz? j + 2x2y72k. 


‘The maximum rate of change is |grad | in the direction of 
grad 


lgrad |” 
lgrad $| = |—2i + 2j + 2k| =2,/3 


At the point (—1, 1,1) we have 


and 
grad 1 Wes ak 
eae a peek 
madd a * al * ak 


(ii) The rate of change of ¢# in the direction of e is (grad p)-e. 


4 
Fore = 31+ $k we have 


(grad p)-e = ee + ane 
5 5 
At the point (2,1, —1) we have (grad ¢)-e = —4. 
5. For two scalar fields @ and y we would expect that 


Vigh) = (Vow + d(VY). 


To check this result we consider V(y) in terms of the 
partial derivatives. 


é 
Foo) = (On + SOW) + Lcd 


é oy 0p op). 
Sy oli (Sy + 0 j 


ow i ba oy i) 


= (VOW + (VY). 


Solutions to the exercises in Section 3 
1 t=0, (a,0,0); 

t=n/4, (a/,/2,/,/2,0); 

t=n/2 (0,a,0); 

t=3n/4, — (-a/,/2,a/,/2,0); 

t=n, (—a,0,0); 

t=3n/2,  (0,-a,0). 


The curve lies in the plane 2 = 0, and is a circle of radius a 
with centre at the origin, 


2. ()x=acos, yy =asine; 


4 =0, t, = 2n; 


w= t(acosti + asintj); 
a 


dr ioe 
qo asl + acos tj; 


2m 
$ wdr= [ Odr=0. 
G 0 


(ii) The parametrization is the same, so now 
u = —asinti+ acostj; 
dr 


—- = —asinti + acos tj; 
a ti+ ‘J 


t, =0, ty = ly 

F = 2ti+ (4 -2)j +k; 
aes 2p. 
ge it 2th + 37k; 


1 
F-dr= f (5t5 — 20dte = —2. 
ce 0 

ae 


4 x= y=?, 


th=1, =0; 

F = 2ti + (t*-2)j+ ek; 
de 

=i 3k; 

a i+ 2rj + 3°k 


fra for —2t)dt = +2. 


5. C, is defined in Exercise 3 


z 


x 


(i) Line integral along C,: 


xem yst rat; 


h=0 h=h 
F=2i+(P —2)j+0k; 


Cre r 
zz i+j+2k; 
1 
f F-dr= f (2+ 30° — 2)dt = 4, 
lez O 
(ii) For the curve C, we go along C, and return along C2, so 


f Fede = fre = [re =a =a: 


3 

6. org=2 4% 

then grad = 201 4 5 4 Ek 
=x i+yj 


and this is the vector field used in Example 4. 
(ii) From Example 4, 


Ra Norge Wea Ca Sera 
Peseta 
= (X25 ¥2522) — POX, V1 21)- 


7. Suppose that A and B are any two points in the domain 
of a scalar field ¢ and let C be any curve joining A to B. 
z, 


B 
Cc 
y 
x A 
Let a parametrization of the curve be 
x=),  y=alt, z= Al), 
with ¢ =f, and =f, at A and B respectively. 


Then grad = 541 + $F) + Eo: 


dr dx. dy, dz 
a al + a + a” 
and hence 
de _apdx | dgpdy | dpdz 
Babb T= Sede * Oydt * deat” 
Using the chain rule from Unit 25, this can be written as “. 


thus 
D a 
Js o-dr = i (eras oS) at 


- [fa 


= b(t) — Oh) 
= $(B) — (A). 


This result is true for all paths joining A and B. Hence, the 
value of the line integral of grad @ depends only on the end 
points of the path (ie, A and B) and not on the path itself. 


Solutions to the exercises in Section 4 


1. (i) The velocity vector at the point P is given in terms of 
i andj by u = —qasin@i + wacosOj 


and acos@ = x 
asind=y 


so that u= —qyi + wxj. 
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(i) x = boost, y = bsint; 
t, =0, t, = 27; 
u= —obsinti + wbcostj; 
dr 


ro —bsinti + beostj. 


Cc 


§ ude = fis (cob? sin? t + wb? cos? t)dt 
c 0 
a8 
-f abdt 
0 
= 2nwb*. 
é 6 
2 s+ 2)-1+ Fr 0+ 
dy)? £"0) 
tee 
y' 5 ay?" 
ant s{y- 2) = s0y-Zr0y+ (2) 
of") |, 
es 


Av 4) -s[v4 9 i 


eee 


Hence 


-f'v) 


_ oF £0) 
4° 6 


and lim 
ay0 


3. (i) For the square given in Figure 3, 


$ war f wdr+ [ wart [ war + [ w-dr. 
f ie ne he pa 


Along AB, r= xi—=j since y= 3 


so taking x as a parameter along here gives 


abe eae eee agi 
urdi =f (r-# d =| ie 
i O* Necays\) (Ge) aleuat2 zi 
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b 
Along BC, ratityi since x= and y will be a 


parameter along here, giving 


2 de "2 wb _ wb? 
-ar= ff (ra =[ ha 
fe Pea Wie ee 8 


sb? 
ude =O 


cob? 
and f ade 
iB 2 Me 2 


Henoe udr = 2wb?, Now area A = b? for the square, so 
é 


Similarly 


ee! ee | 
lim fw de = lim 5 (2wb*) = 2u. 


ano 
f ude + [ wde+ [ wde+ [ u-dr. 
4B Be cp Da 


Gi wd 


Along AB, r = xi+ (rs - ji so 


J ude = ~ofr — 3) ax 
AB xn = W/2 2 
b| 
oaliedl 


Along BC, r= (x + °) i+ yj, so 


yo b/2 


f ude = afro +5) = o6(x0-+ 4), 
ac yo b/2 2, 2 


Along CD, r= xi-+ (49 +3), 50 


ova b ’ 
fowe= fl fro +§} 4x = (ro +4). 
co Perr 


Along DA, r = (x - ;)i + yj, so 


[ee [efor t 


Thus f u-dr = 208 so that tim +4, u-de = 20 again, 
fi anode 


4. (i) ey = —sin Oi + cos0j 
y 


x 


and sin@ = 
Ir} 


and cos@ = 


Vie Se 
H Sey ee EE 
aaa oh eee are 


47 


(i) u= fe so using the result in part (i) we obtain 
u - 4-3 - *i) (where d = /x? + 7) 


¢ a 
=pl-yit+ xj). 


(iii) Suppose that AD and BC are arcs of circles of radii a 
and b respectively. 


y 


fwde= f wide + f ude + f wde+ f usdr 
f an be ce ns 


limsn 
Fe 


Similarly for the side CD, u-dr=0. 


cD 
Along BC, we can use the parametrization 
x=bcost, y=bsint 


andt,=¢,%=¢+a. 


On substitution w = “S(—bsinti + beosj), 


dr 4 
a —bsinti + becosej, 
dr alt IR 
so that wane and f (e flar= ce 
Similarly, along DA, 
x=acost, y=asint 
and (; = + a, [: = ¢; so 


f urdr = edt = —ca, 
DA J 


ot 


Hence u-dr = 0+ cx +0 ~ca=0. 
é 


Gy) (curly), = tim + usd =O since u-dr= 0. 


a0 Ac 
5. (curlw), = -F 


@ = —ayi+ oxj 


so—-= — = 


ax ay 
so (curl u). = 2m, 


’ -c 
Oia Fim 
ye ext = x 
os G+) Carr Ory’ 
oF =e 5 ey? _W=x 
y +) +P &+y¥P 
so (curlu), = 0, 
6 (@)For=xit+yj+zk; 
R=x hay Kan 
curlF = 0. 
(i) F = fre; 
=Si0=,  R=fO  Fr=S= 
i component of curlF is 
Oa _ OF 4 (file _ (ena 
dy dz dr\ r dy ar\ r Jaz 


(using the Chain Rule) 


L\? <0, 


Similarly the j and k components of curlF are zero, so 


curl (f(r)#) = 0 for all functions f of r. 
(iii) F = (y? + 22) + (xy + 62)j + (2? + 2x2 + y)ks 
oF, oF ow 
i component of curlF is ee =1-6= 
j component of curlF is oF 8 9 22 = 2-2); 
dz Ox. 
OF, 0, 
k component of curl F is = - ES =y-ly=-y; 
so curlF = —Si + 2(1 — z)j — yk. 
(iv) F = grad = 148 + e+ e a: ki 
ar, oF _ a (26 aa 
dy dz ady\az] az “3 i 
oa _e od) 08 oh 
dz ox dz\dx} ax\az} 
a _9F _ 0 (26) _ 0/06) _ 9, 
ox ay a ay] ay\ox}” 
hence curl (grad ¢) = 


1. VxF= (2 widens] x (Fit Fij + Fk) 
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and using the formula for the cross product given in Unit 14 
this becomes 


[of _ 2) , (2h 2h) ea 
ay ae) Vee ox ax ey 
=anlF. 


8. (i) F 
curlF 


2xi + (xz — 2)j + xy kj 
— x) + (0 — y) + k(z — 0) 


= —yj+ zk. 
Since curlF # 0 in general, F is non-conservative. 


(ii) From Exercise 6(ii) curl( (rt) = 0 for all functions f(r). 
For G = r"t, f(r) =r" so that curl G = 0 and hence G is 
conservative 


9, If F is conservative then F-dr =0 for all closed curves 


ts 
C. Consider a closed curve C passing through the two points 
Aand B. 


$ Fede = f Fede+ f F-dr=0. 
td APB BQA 


Thus f Fede= —f Fedr. 
APD BQA 


But f Fedr= — Fedr; 
BOA age 

Fedr= 

Are AQ 


hence F-dr and this is true for any path 


between A and B. So the line integral ) Fedr is 


AB 
independent of the path, and depends only on the end points 
Aand B. 


Cott 


